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1. Motivation 

Recently the A-parameter of two-flavour QCD has been computed by the ALPHA collabo- 
ration in the Schrodinger functional scheme in units of a low-energy scale L max [jl]]. The latter 
scale is defined implicitly by the renormalized coupling in that scheme taking a particular value, 
f| F (/i = 1/L) = 4.61. The determined quantity AsF^max is a continuum, universal result. The 
connection between Asf and A^j§- being known exactly [||], what remains to be done, for the re- 
sult to be phenomenologically useful, is to trade L max for an experimentally accessible low-energy 
scale. This was done provisionnally §\§ using the chirally extrapolated Sommer reference scale 
i"0 [|3^ [ljl, however choosing a more directly accessible physical quantity will reduce the systematic 
uncertainty. Our present goal is to use the kaon decay constant Fk as low-energy scale. This will 
also allow us to reduce the overall uncertainty on our recent determination of the strange quark 
mass ft 

More specifically, a possible strategy is to compute L max F ps at a tuned quark mass such that 
^ (m re f , m re f ) = TjHexp- Neglecting the quenching of a third quark of mass m re f, SU(3) chiral 
perturbation theory then connects F ps (m Te f,m ie f,m ie f) to F ps {m u ,mii,m s ), which is equated to its 
experimental value of 113MeV. At one-loop level [Q], with /ip = Jfl F i log{Mp/}X 2 } and Fq = 
88MeV we have 

{m +Z s )/2 = 1 + l^-TMr+^ + f Mn)-L 6 (n)].[M 2 K -Ml] (1.1) 
= 1 + 0.10 + 0.47-10 3 [L 4 (M p )-L 6 (M p )], 
F p s (m M ,m ref ) = 1 + | (At7t _ 2 ^ + ) + ^. L4(At) . [M 2_ M 2 ] (12) 

= 1 + 0.051 + 0.12-10 3 L 4 (M p ), 

where phenomenological studies point to 10 3 |L 4 | < 1 and 0.2 < 10 3 [L 4 (M p ) -L 6 (M p )} < 0.4 
These ranges can be narrowed down by lattice simulations. Using the tree-level relation \M\ — 



= M\, it can be verified that the dependence on \i cancels in Eq. (LI) and (|L 
Although the original purpose of the Schrodinger functional was to provide the setting for 
the implementation of a finite-volume renormalization scheme, it was shown [§§] in the quenched 
approximation that it also provides a competitive tool for spectrum and matrix element calculations. 
Roughly speaking, the Dirichlet boundary conditions in the time directions are now exploited in 
the same way as in the source method that was used in the early days of glueball spectroscopy [Hj. 



2. Simulations at m re f 

We use the Aff = 2 O(a) improved Wilson formulation with plaquette gauge action S g . Our 



simulations are performed with the Hybrid Monte-Carlo algorithm with (symmetric) even-odd [ 1C] 



and Hasenbusch preconditioning [pi]]. The action thus reads S = S p f + Sdet + Sg> with 
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Figure 1: Comparison of the pseudofermion and determinant forces in the time-slice T/2 with symmetric 
and asymmetric even-odd preconditioning on an 8 4 lattice at j3 = 6.32 (n = 1 here, and the normalization of 
the force is as in [ 18 1). 



where the pseudofermions (PFs) 'live' on the odd sites only and, in the notation of [10], 

Qa = 7s (M 00 - M oe M^M eo ), Q = M- a l Q A , (2.2) 
Sdet = (-2)logdetM ee + (-2)logdetM 00 , (2.3) 

The parameters (Ok,Pk) satisfy 

°k=Pk-i-Pk, k = l,...n-l; 0,-1=0. (2.4) 

The Hermitian matrix Qa is thus the Schur complement of Q = y$D with respect to even-odd pre- 
conditioning. Note that Q on the other hand is not Hermitian. The matrices M DO ,M ee are diagonal 
with respect to space indices, and their determinants are thus represented exactly on each configu- 
ration. All the simulations presented here are with n = 2 unless otherwise stated. 



2.1 Recent algorithmic improvements 

The advantage of the symmetric form of even-odd preconditioning was first demonstrated 
in [jl2|] with n = 1. See [|l~3|] for a study with n = 2 (where the shift p was however implemented 
on D rather than Q). We find that on a 16 4 lattice at j8 = 5.2, the stepsize 8z can be increased 
from 1/28 to 1/20 at a constant acceptance above 90%, when switching from the asymmetric to 
the symmetric form. This is easily explained by a reduction in the magnitude of the force, as 
illustrated on Fig. |1J for a small system. The operator Q is better behaved in the ultraviolet than Qa ; 
the improvement is thus not expected to depend strongly on the size of the system. Secondly, in 
the Hamiltonian and force computation, the equation (Q A + p 2 )y = (j) arises, to be solved for \y. 
We then proceed as follows, with M^, stored: 

0=iH^V, (Q-ipM- 1 )(Q f + ipM- 1 )ijf = $, y = M- l y. (2.5) 
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Table 1: The figure-of-merit v[0] for the plaquette, the effective pseudoscalar mass and decay constant at 
j3 = 5.3 in the range of masses m s /2 < m < m s on 24 3 x 32 lattices. 



We use the conjugate gradient algorithm for the inversion involving Q, and on the 16 4 lattice find 
about a 20% reduction in the number of iterations necessary to solve the equation with p = for 
a given target accuracy with respect to solving the equation involving <2/t- Thus the symmetric 
variant is also a superior solver preconditioner as compared to the asymmetric one. 

It was recently suggested to introduce a separate step-size 8tk for each pseudofermion [14, 15] 
(a smaller step-size for the gauge force has been in use for a long time [16]). With n = 2 one then 
has one additional parameter to tune, and it turns out that the simultaneous optimization of po and 
5to/5ti leads to significantly smaller values of po j fOj [l7| ] than if one constrains 5 To = 5X\. We 
find a gain resulting from the introduction of the extra parameter of about 1.5 on a 24 3 x 32 lattice 
at m ~ m re f, while no observable improvement could be obtained on the 16 4 , j8 = 5.2 lattice. This 
suggests that the gain increases with the condition number of Q. 

We have also explored 'hybrid' integrators, where the leap-frog integration scheme is used for 
the force associated with 0i and the Sexton-Weingarten [ ff6| ] for that associated with <po, and found 
the performance to be comparable. The motivation is to benefit from the robustness of the leap-frog 
integrator for the most irregular force to insure stability, while using a higher order scheme for the 
smoother forces. 



3. Trajectory length dependence of autocorrelation times 



We report some results on the algorithm figure-of-merit v introduced in [18], which we gene 



ralize to any observable O. If Tj nt [O] is its autocorrelation time in units of trajectories, we have 

V [O] = KT 3 • Aff£ lver calls • Tin, [O] . (3.1) 

By 'effective number of solver calls', we mean the number of times the equation Q\y = has to 
be solved per trajectory, plus the number of times the equation {Q 2 A + pfyy = has to be solved, 
weighted by the number of solver iterations this takes relative to the former equation. It turns out 
that these two contributions are roughly equal. The values of v for the plaquette, the effective 
pion mass and pseudoscalar decay constant in the middle time-slice are given in Tab. |I[ There 
is certainly no sign that the algorithm's performance worsens as the quark mass decreases in the 



explored region. The figures for the plaquette P are very similar to the ones quoted in [19], and 
also in [ fl5] ] (if one uses the same definition of A^g' ver calls ). 
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Figure 2: Autocorrelation function of the correlators f\ and f P (Nf = 2, j3 = 5.3, K = 0.1355, V — 24 3 x 32), 
for different trajectory lengths. 

At K = 0.1355, two different runs with trajectory length 1/2 and 2 were performed. While 
the plaquette actually 'prefers' the shorter trajectory, the observables of interest (especially / ps ) 
have a smaller autocorrelation time with the longer trajectory. Determining autocorrelation times 
accurately is difficult, and therefore it is useful to look directly at the normalized autocorrelation 
function p(t). In Fig. || this function is shown for two fermionic correlators: f\ which corresponds 
to the propagation of the (ud) quarks from one boundary to the other and fp, where the pseudoscalar 
density in the middle time slice annihilates the states created by the boundary (see [Jsj] for the precise 
definition of these correlators). A noticeable difference is seen at small t. To make the observation 
more significant, we repeat the comparison in the quenched theory with an equivalent quark mass, 
using the HMC algorithm for the pure gauge update, and with a much smaller spatial volume (see 
Fig. |H). Here the statistics is high enough to allow us to see a clear difference between the choices 
T = 1/2,2 and 4 for the trajectory length. The longer trajectories are undoubtedly superior in this 
situation. 

The reduced autocorrelation times observed with long trajectories should be reflected in faster 
thermalization of the same observables. We thus perform the following experiment. We start from 
an ensemble of 16 small- volume and small-quark mass configurations with a background chromo- 
electric field. The correlators fp and f' p , defined by the correlation from one or the other boundary 
to the middle time-slice, are thus not equal. We then change the SF boundary conditions so as to 
remove the background field, after which fp = f P holds exactly in the new thermalized ensemble. 
Fig. |] shows how fast f P — f' p thermalizes to zero with two different choices of trajectory length. 
The advantage of z = 2 over % = 1/2 is obvious in the early stages of thermalization. 

4. Conclusion 

We are performing large-volume Nf = 2 simulations with 0(a) improved Wilson quarks in 
the Schrodinger functional at m ~ m s /2 to set the low-energy scale for our determination of the A 
parameter [jl|] and the strange quark mass [Qj. Apart from the L/a = 24 simulations presented here, 
simulations at L/a = 32 and larger j8 are underway. 



5 



Trajectory Length and Autocorrelation Times 



Harvey Meyer 



Normalized autocorrelation of f1 



+ T=1/2 

o x=2 

T=4 



10 20 30 40 50 60 70 80 90 100 

MD time 



Normalized autocorrelation of fP 



0-8 I 

m It 

o.6- a 



+ t=1/2 
o x =2 
x=4 



1¥Pt 



10 20 30 40 50 60 70 80 90 100 

MD time 



Figure 3: Autocorrelation function of the correlators f\ and fp in quenched QCD (j3 = 6.0, K = 0.1338, 
V = 8 3 x 32), for different trajectory lengths. 
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Figure 4: Thermalization of fp — fp, whose expectation value vanishes. The data points are averaged over 
16 configurations taken from an ensemble where the expectation value is non-zero. 



We have shown that the trajectory length z is an important parameter of the Hybrid Monte- 
Carlo algorithm. In practice we have found that z > 2 is never a worse choice for long-distance 
observables based on fermionic correlators. The reduction in autocorrelation time is in some cases 
as large as a factor two. We have also provided evidence that slow-thermalization problems can be 
alleviated by choosing longer trajectories. For more details, as well as a discussion of stability and 
reversibility issues, we refer the reader to [SJ. 
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